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Partial dynamical symmetry describes a situation in which some eigenstates 
have a symmetry which the quantum Hamiltonian does not share. This prop- 
erty is shown to have a classical analogue in which some tori in phase space are 
associated with a symmetry which the classical Hamiltonian does not share. 
A local analysis in the vicinity of these special tori reveals a neighbourhood 
of phase space foliated by tori. This clarifies the suppression of classical chaos 
associated with partial dynamical symmetry. The results are used to divide 
the states of a mixed system into "chaotic" and "regular" classes. 
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Symmetry plays a central role in affecting the character of classical and quantum dynam- 
ics. Constants of motion associated with a symmetry govern the integrability of the system. 
Often the symmetry in question is not obeyed uniformly, i.e. only a subset of quantum 
states fulfill the symmetry while other states do not. In such circumstances, referred to as 
partial symmetries, the symmetry of these "special" (at times solvable) eigenstates does not 
arise from invariance properties of the Hamiltonian. Selected examples in this category are 
adiabatic regular states in the stadium billiard [Q; regular quasi-Landau resonances of a hy- 
drogen atom in strong magnetic fields |§; discrete nuclear states embedded in a continuum 
of decay channels ||; partial SU(3) symmetry found in deformed nuclei ||]. Hamiltonians 
with partial symmetries are not completely integrable hence can exhibit stochastic behavior. 
As such they are relevant to the study of mixed systems with coexisting regularity and chaos 
||, which are the most generic. 

The effect of discrete symmetries on a mixed-phase-space system was examined in ||, 
while the consequences of breaking a discrete symmetry on the spectral statistics of a chaotic 
Hamiltonian were discussed in [[7|. In the present paper we focus on continuous symmetries 
(associated with Lie groups) which can be conveniently studied in the framework of alge- 
braic models. Such symmetry-based models are amenable to both quantum and classical 
treatments and have been used extensively in nuclear and molecular physics |§ . Their inte- 
grable limits are associated with dynamical symmetries in which the Hamiltonian is written 
in terms of Casimir operators of a chain of nested algebras. The eigenstates and eigen- 
values are labeled by the irreducible representations (irreps) of the algebras in the chain, 
and are known analytically. Both the quantum and classical Hamiltonians are completely 
solvable ||. Non-integrability is obtained by breaking the dynamical symmetry and may 
lead to chaotic dynamics ||[10|. To address situations of mixed dynamics, an algorithm was 
developed JET]] for constructing algebraic Hamiltonians with partial dynamical symmetries. 
Such Hamiltonians are not invariant under a symmetry group and yet possess a subset of 
"special" solvable states which do respect the symmetry. In the context of a nuclear physics 
model involving five quadrupole degrees of freedom, it was shown that partial dynamical 
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symmetry (pds) induced a strong suppression of classical chaos [ I2[ . This was true even 
though the fraction of special states vanished as h 2 , so one might have expected no classical 
effect. In order to better understand this effect, we consider a simpler model and use its 
pds to infer relationships between the classical and quantum dynamics of a Hamiltonian in 
a mixed KAM regime. 

As a simple test-case we consider a model based on a U (3) algebra. This algebra (with 
fermionic operators) was considered previously in the context of chaos, but not regarding 
pds Jul- Here we employ a realization in terms of three types of bosons , b\ satisfying 



the usual commutation relations. The nine number- conserving bilinear products of creation 
and destruction operators comprise the U(3) algebra. The conservation of the total boson- 
number = h a + n b +h c [h a = a^a with eigenvalue n a etc.) ensures that the model describes 
a system with only two independent degrees of freedom. All states of the model are assigned 
to the totally symmetric representation [N] of £7(3). One of the dynamical symmetries of 
the model is associated with the following chain of algebras 

U(3) d 17(2) D 17(1) (1) 

Here £7(2) = £77(2) x U a b(l) with a linear Casimir n ab = h a + n b (which is also the generator 
of U ab (l) ). The generators of 577(2) are J + = tfa, J_ = alb, J z = (n b — h a )/2 and its 
Casimir J 2 = n ab (h ab + 2)/4. The subalgebra U(l) in Eq. ([]]) is composed of the operator 
J z . A choice of Hamiltonian with a £7(2) dynamical symmetry is 

H = uo a a)a + u b b^b = n ab — 2AJ Z (2) 

where uo a ^ b — 1 ± A, and A is introduced to break degeneracies. Diagonalization of this 
Hamiltonian is trivial and leads to eigenenergies E na)Tlb = u a n a + ujbn b and eigenstates 
\n a , n b , n c ) or equivalently \N,J,J Z ) where the label J = n ab /2 identifies the SU{2) irrep. 
These are states with well defined n a , n b and n c = N — n a — n b . To create a pds we add the 
term 

Hi = b^(b^a + b ] c + a)b + c ] b)b , (3) 



which preserves the total boson number but not the individual boson numbers, so it breaks 
the dynamical symmetry. However states of the form \n a , = 0, n c ) (or equivalently 
\N, J = n a /2,J z = —J) ) with n a = 0, 1,2, ...N are annihilated by Hi and therefore 
remain eigenstates of Hq + BH\. The latter Hamiltonian is not an SU(2) scalar yet has a 
subset of (N + 1) "special" solvable states with SU{2) symmetry, and therefore has pds. 
There is one special state per SU{2) irrep J = n a /2 (the lowest weight state in each case) 
with energy uo a n a independent of the parameter B. Other eigenstates are mixed. Although 
H and Hi do not commute, when acting on the "special" states they satisfy 

H , Hi\\n a ,n b = 0,n c ) = 0. (4) 
To break the pds we introduce a third interaction 

H 2 = a f c + c f a + 6 f c + c f 6 . (5) 

The complete Hamiltonian is then 

H = H + BHi + CH 2 . (6) 

For B = C = we have the full dynamical symmetry; for B ^ 0, C = we have partial 
dynamical symmetry and for C ^ we have neither. 

The classical Hamiltonian 7i c i is obtained from (§) by replacing (a^", 6^, c^") by complex 
c-numbers (a*, (3*, 7*) and taking iV — > 00. The latter limit is obtained P,|T0|| by rescaling 
B = NB, a — > a/y/N etc. and considering the classical Hamiltonian per boson Ti. = 7i c i/N . 
In the present model the latter has the form 

H = H + BHi + CH 2 ■ (7) 

Number conservation imposes a constraint a* a + (3* (3 + 7*7 = 1, so that the phase space is 
compact and four-dimensional with a volume 2n 2 . The total number of quantum states is 
(N + l)(iV + 2)/2. Assigning, to leading order in N, one state per (2nh) 2 volume of phase 
space, we identify % = 1/N, so that the classical limit is iV — > 00. 



In all calculations reported below we take A = 0.8642 and N = 60. As a first step, we 
fix B = 0.5 and vary C. For the classical analysis we randomly sample the phase space 
and determine the fraction fi of chaotic volume. For the quantum analysis we evaluated 
the energy levels, calculated the nearest neighbors level spacing distribution of the unfolded 
spectrum and fitted it to a Brody distribution |0 . The Brody fit parameter u> is expected to 



be for integrable systems (Poisson) and 1 for chaotic systems (GOE). As shown in Fig. 1, 
both of these measures indicate a suppression of chaos near C = similar to the results of 
Ref. Jl^|. To appreciate the strong effect of the pds (at C = 0) on the underlying dynamics, 
it should be noted that the fraction of the solvable states \n a , rib = 0, n c ) is 2/ (N + 1), which 
approaches zero in the classical limit. To measure the extent to which each eigenstate 
has SU(2) symmetry, we define variances of = ( 1 J/|n?| 1 I r ) — (\l/|nj|\l/) 2 (i = a,b). A state 
which belongs to just one irrep of SU(2) (with well defined J, J z ) has zero variances, while 
a mixed state has large variances. These variances have the same physical content as the 
entropies considered in Ref. ||12|| . It is instructive to display the average (n a ) and variance 
of each state, as done in Fig. 2. SU{2) pds is present in Fig. 2a (B ^ 0, C = 0), Fig. 2b is a 
blow up of Fig 2a and in Fig. 2c the symmetry is completely broken (C ^ 0). In Figs. 2a-b 
we see states with zero variance. These are just the special N + 1 states (rib = 0) discussed 
before, which preserve the SU(2) symmetry. In addition, we see families of states with small 
variance and small (rib) which suggests that the presence of partial symmetry increases the 
purity of states other than the special ones. By contrast, in Fig. 2c we see no particular 
structure because of the destruction of the pds for C ^ 0. 

Considerable insight is gained by examining the classical phase space structure in terms of 



action-angle variables a = \fj~ a exp(— i6 a ), (3 = \fJb~ exp(—i6 b ) and 7 = \fT c = y/1 — J a — J a . 
The 9 a = —tt/2 Poincare section is shown in Fig. 3 for energy E — 1.0. When SU(2) pds is 
present (B ^ 0, C = 0) we see in Figs. 3a-b a torus with J b = 0, and additional perturbed 
tori in its neighborhood (small Jb). This structure is absent when the symmetry is completely 
broken (C^0), as shown in Fig. 3c. The features in Fig 3 persist also at other energies. To 
understand them, we recall that for B = C = 0, the Hamiltonian (|7|) is integrable and all 



trajectories wind around invariant tori. By standard torus quantization (without turning 
points) the actions are quantized as Jj = rtfi = Ui/N (i = a,b). In the integrable limit 
quantum states are associated with toroidal manifolds in phase space. In case of a partial 
symmetry (B ^ 0, C = 0) we are led by analogy with Eq. (^) to seek manifolds M. in phase 
space on which 

{Wo,Wi}| A4 =0 (8) 

vanishes even though the Poisson bracket is not zero everywhere. In addition, we demand 
{{7io , 7~Li}, 7~(-o+7~(-i}\m — (in analogy to the quantum relation [[ H , Hi] , H +Hi] \n a , n b = 
0,n c ) = 0) so that a trajectory starting on M. remains on M.. The solution to these 
conditions is the manifold J b = (3* (3 = 0, which may be interpreted as a (degenerate) torus 
of the Hamiltonian. It is also a stable isolated periodic orbit of Tio + BTCx. Quantization 
of the torus with J b = proceeds exactly as before, so we correctly predict no change in 
the quantum energies associated with it. The manifold M. (J b = 0) is the direct classical 
analogue of the special quantum states \n b = 0). It refers, however, to a region of phase 
space of measure zero, and so cannot by itself explain the observed (global) suppression of 
chaos. However, as suggested by Fig. 3, the presence of pds induces a quasi-regular region 
foliated by tori in the vicinity of the special torus. We can understand the dynamics on a 
finite measure of phase space by performing a perturbative calculation in the neighbourhood 
of M. 

For the classical perturbation calculation we set C = in Eq. (0) and treat B as an 
expansion parameter, assuming BTii in Eq. (0) is small in the neighbourhood of the special 
periodic orbit. We seek a generating function for the canonical transformation (J, 6) — > 
(I, (f>) such that the Hamiltonian depends only on the new actions I and not on the new 
angles cf>. The first-order result in standard classical perturbation theory is 

S{I, 6) =16- 2Bll' 2 



iy 2 sm(9 a -e b ) l¥ 2 sin9 b 



(Auj) u b J ' (9) 

where Atu = uj a — uj b and J = |jj , cf) = . There is no first order shift in the energy hence 
Tt(I) = oo a I a + oo b I b and the calculation is valid for BI b ^ 1. The second order correction 



reproduces well the perturbed tori on the Poincare sections as shown in Fig. 3b. The classical 
variances are defined by averaging with respect to the angles (9 a , 6b), e.g. a 2 d = (J%) — ( J a ) 2 - 
To the same order in perturbation theory, (J a ) = I a , (J b ) = I b and the variances in J a and 
Jb are 
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The same calculation can also be done in quantum perturbation theory with the result 
B 2 n b 



(10) 
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(Alu) 2 (tub) 2 ) y " ' \(Alu) 2 (cub) 2 
In the semiclassical limit, n a , rib, n c ^> 1, Eqs. (|10"D and (|TTD agree to leading order. In 
Fig. 2b we show the results (denoted by (+) ) of the quantum perturbation theory (to order 
B 5 ). We see that the first few families of states are reproduced. It is these states which we 
can recover from perturbation theory and whose approximate symmetry is induced by the 
symmetry of the special states. 

The following physical picture emerges from the foregoing analysis. Near the special 
orbit, there are KAM tori, some of which are quantized. The quantum eigenstates lie on 
these tori, so knowing the classical variance of the actions of the tori tells us the variances 
of the states themselves, in the semiclassical limit. Large variances indicate the extent to 
which the corresponding states fail to respect the symmetry. This provides a measure for a 
separation of regular and irregular levels, as conceived by Percival |L4j] and which recently 
gained numerical support |6|JT5||. In the present model the quantum states can be grouped 
into three classes: i) the special states, which observe the symmetry; ii) the "almost special 
states" which are accessible by perturbation theory; iii) the rest of the states, which are 
mixed. As in 0, the frontier between regular states (sets (i) and (ii) ) and irregular states 
(set (iii) ) is not sharp. 

In summary, we have considered the effect of partial dynamical symmetry (pds) on the 
quantum and classical dynamics of a mixed system. At the quantum level, pds by definition 



implies the existence of a "special" subset of states, which observe the symmetry. The pels 
affects the purity of other states in the system; in particular, neighboring states, accessible 
by perturbation theory, possess approximately good symmetry. Analogously, at the classical 
level, the region of phase space near the "special" torus also has toroidal structure. As a 
consequence of having pds, a finite region of phase space is regular and a finite fraction of 
states is approximately "special". This clarifies the observed suppression of chaos. Based 
on these arguments and the numerical results of Ref . (T^j , we anticipate the suppression of 
chaos to persist in higher dimensional systems with partial dynamical symmetries. 
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FIGURES 

FIG. 1. Classical (//) and quantum (oo) measures of chaos (denoted by (•) and (x) respectively) 
versus C for the Hamiltonian (6) with B = 0.5. 

FIG. 2. The values of (n a ) and of the variance o\> (denoted by o) of each eigenstate of the 
Hamiltonian (6). (a) B = 0.5, C = (partial dynamical symmetry), (b) a blow up of (a) with 
superimposed results (denoted by (+) ) of quantum perturbation theory. The families of states 
with low o i) have small values of {rib). ( c ) B = 0.3, C = 0.5 (broken symmetry). 

FIG. 3. Poincare sections Jb versus 9b at energy E = 1.0. (a) B = 0.5, C = (partial 
dynamical symmetry), (b) a blow up of (a) with superimposed results (dashed curves) of classical 
perturbation theory, (c) B = 0.3, C = 0.5 (broken symmetry). 
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